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The hybrid relay selection (HRS) scheme, which adaptively chooses amplify-and-forward (AF) 
Q , and decode-and-forward (DF) protocols, is very effective to achieve robust performance in wireless 



networks. This paper analyzes the frame error rate (FER) of the HRS scheme in general cooperative 
wireless networks without and with utilizing error control coding at the source node. We first develop 
an improved signal-to-noise ratio (SNR) threshold-based FER approximation model. Then, we derive 
an analytical average FER expression as well as an asymptotic expression at high SNR for the HRS 
O . scheme and generalize to other relaying schemes. Simulation results are in excellent agreement with 

the theoretical analysis, which validates the derived FER expressions. 
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Index Terms 

Hybrid relay selection, frame error rate, SNR threshold-based approximation, cooperative com- 
munications. 

I. Introduction 

Cooperative relaying has been shown to be an effective technique to improve the system 
performance in wireless networks by allowing users to cooperate with each other in their 
transmissions. Two of the most typical cooperative schemes are amplify-and-forward (AF) 
and decode-and-forward (DF). The performance of AF protocol is mainly limited by the noise 
amplified at the relay during the forwarding process, especially at low SNR. The performance 
of DF will be degraded when the relay fails to decode the received signals correctly and the 
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process of decoding and re-encoding will cause serious error propagation. To overcome these 
limitations, various cooperative schemes have been reported, such as signal-to-noise-ratio 
(SNR) threshold based selective DF [[Il-Hl, cooperative-maximum-ratio-combining based 
DF [El, ja, decode-amplify-forward (DAF) [Hi, [H, link adaptive relaying DF M-W and 
log-likelihood-ratio (LLR) threshold based selective DF [[HI, [[Bl, etc. 

Recently, the hybrid relaying protocol (HRP) has received a lot of attention. It adaptively 
combines the merits of both DF and AF by forwarding clean packets in DF if decodes 
correctly and forwarding soft represented packets in AF if decodes incorrectly) [[T4ll - 
ifTSl . Its performance can be further improved by incorporating relay selection in HRP, 
which is referred to as hybrid relay selection (HRS) [IT4ll - [l24ll . In the HRS scheme, for 
each transmission, all relays are divided into two groups, a DF group and an AF group. A 
relay is included either into the DF group if it decodes correctly or into the AF group if it 
decodes incorrectly. The destination node selects a single optimal relay node with maximum 
destination SNR from either the DF group or the AF group to forward the packet. 

Although the principle of the HRS scheme is simple, the calculation of the analytical 
FER is non-trivial even for single-relay cooperative networks. Therefore, most earlier work 
have been focusing on the approximate average FER analysis. Nonetheless, still very limited 
results have been reported in the literature so far. In [16|, Souryal analyzed the FER of 
single-relay networks using hybrid forward scheme in block Rayleigh fading channel. In 
IfTTlj . |[2T]| . Li proposed the hybrid relay selection scheme by combining hybrid forwarding 
and relay selection. In [[22ll . Song extended the HRS scheme with differential modulation. 
Only very loose upper bounds were provided for the FER of the HRS scheme. In [25], the 
authors mainly analyze the FER of the all-participate scheme without relay selection, where 
all relays participate in forwarding signals from the source. 

In this paper, we derive the analytical FER of the HRS scheme in general cooperative 
wireless networks with and without applying convolutional coding at the source node. The 
contribution of this paper is threefold: 1) we develop an improved SNR threshold-based 
FER approximation model, which is simple and accurate for general diversity systems, and 
obtain the SNR threshold in an analytical expression; 2) we derive an analytical approximate 
average FER expression, and its simplified asymptotic FER expression at high SNRs for the 
HRS scheme; 3) we generalized the derived theoretical FER expressions to other relaying 
schemes, e.g. the AF-RS scheme and the PDF-RS scheme. 

The rest of the paper is organized as follows: We describe the system model and the HRS 
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scheme in Section [III In Section [nil we develop an improved SNR-threshold based FER 
approximation model based on cumulative distribution functions (CDF). In Section |IVl we 
derive the analytical approximate and asymptotic FER expression at high SNRs of the HRS 
scheme. In Section |Vl we present some numerical simulation results. And in Section |VIl we 
draw the main conclusions. 

Notation: Boldface lower-case letters denote vectors. Zn represents the n-dimensional 
binary space {0, 1}". For a random variable X, Pr{-) denotes its probability, fx{ ) denotes 
its probability density function (PDF), Fx{-) is its CDF, and E[X] represents its expectation. 
X ~ CJ\f{0, ^7) is a circular symmetric complex Gaussian variable with a zero mean and 
variance Q. Q{x) denotes the Q-function, i.e. Q{x) = e~*^/^ dt. 
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Fig. 1. Diagram of the cooperative wireless system with the HRS scheme. 



II. System Model and The HRS Scheme 

In this paper, we consider a general 2-hop relay network, as shown in Fig. [H consisting of 
one source node, n relay nodes and one destination node. The link between any two nodes 
is modeled as a block Rayleigh fading channel with additive white Gaussian noise (AWGN), 
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where the fading coefficients of the channels are fixed within one frame and vary indepen- 
dently from one frame to another. 

ho, hi^i, and /i2,j are the fading coefficients of the channels from the source node to the 
destination node, from the source node to the i-th relay node, and from the i-th relay node 
to the destination node, respectively. We assume ho ~ CN{Q, Vto), hu ~ CN{Q, ^li), and 
/i2,i ~ CjV(0, ^21)- Similarly, no, riu and n2i are the corresponding additive Gaussian noises. 
We assume no ~ CN'(0,No), nu ~ CN'(0,No), and n2i ~ CN'{0,No). Without loss of 
generality, we consider that all the nodes transmit with the same power E. In addition, all 
channel state information (CSI) needed for decoding is available at the relay nodes and the 
destination node. 

We assume that each terminal in the network is equipped with a single antenna working 
in the half -duplex mode. We consider an orthogonal transmission scheme in which only one 
terminal is allowed to transmit at each time slot. Therefore, one frame transmission in the 
HRS scheme consists of two phases: 



A. The First Phase 

The source node broadcasts the transmit signals with length L, denoted by s, to both the 
destination node and all of the relay nodes. The received signals at the destination node and 
the i-th relay node, denoted by yo and y^j, are Jq — y/EhoS + no, and y^j = VEhuS + n^, 
respectively. Then the corresponding instantaneous SNRs can be calculated as 

E(|noP) =iVo'^°' 

and 

g 

iVo 

respectively, where E(|sp) = L, E(|np) = LNo, and 7 is defined as 

-A A 

B. The Second Phase 

In the second phase, every relay node first decodes it's receive signals. Then it is assigned 
to either the AF group Qaf or the DF group Qdf according to its CRC checking result. 
If the CRC checking result is correct, the relay node is included into Qdf, otherwise, it is 
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included into Gaf- As a result, the instantaneous destination SNR of the link through the 
i-th relay node, denoted by can be expressed as 



if i e Qaf 



7.= r"^'"''" (1) 
72i, if i G Gdf , 

where 72^ represents the instantaneous SNR of the link between the i-th relay node and the 

destination node, which is given by 



S 



l2i = Trl^2iP = l\h2i\'^- 

iVo 

Then, the destination node selects the relay node with maximum SNR 7^ to transmit at 
the second phase 

7^ = max{7i}. (2) 

l<j<n 

The selected relay node m transmits through AF if m G Gat using DF if m G 
GvT- Finally, the destination combines the signals received at both phases by maximal ratio 
combining (MRC) to decode. After MRC, the effective SNR of the received signal, denoted 
by 7i^i?5, is calculated as 

'IHRS = 70 + 7m- (3) 

The selection process in the HRS scheme is exactly the same as in the conventional AF 
or DF selection schemes and does not add any extra complexity in system implementation. 
The only requirement for the HRS scheme is that each relay needs send one bit indicator to 
inform the destination if it uses AF or DF protocol. Destination then calculates the overall 
received SNR from each relay accordingly and selects the best relay with largest SNR. 

For simplicity, we define Aq, Aij and \2i as: 

7"o 7"ii 7"2i 

in. An Improved SNR Threshold-Based FER Approximation Model 

In this section, we develop an improved SNR threshold-based FER approximation model 
for general diversity systems. 

We first describe the SNR threshold-based FER approximation model in Subsection IIII-Al 
in Subsection IIII-B[ some earlier approaches are presented, and then we propose an improved 
criterion to calculate the SNR threshold in Subsection IIII-CI 
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A. Introduction of the SNR threshold-based FER approximation model 

The average FER over a block fading channel, denoted by P/, can be computed by 
integrating the instantaneous FER over AWGN channel, represented by Pj^i^j), over the 
fading distribution [|26ll 

/■oo 

Pfi^)= / Pf(7)/7(7,7)d7, (5) 
Jo 

where 7 and 7 denote the instantaneous and average SNR, and /^(■) denotes the PDF of 7. 

Although ([5]) is an exact expression for Pf, its closed-form expression is difficult to 
evaluate. By assuming the instantaneous FER is 1 when instantaneous SNR 7 is below a 
threshold 7^, otherwise it is 0: 

^^^(717 < 7t) ~ 1 and (7I7 > It) ~ 0, 
the average FER can be rewritten as 

Pf= I Pf(7)/,(7,7)d7+ / Pf(7)A(7,7)d7 

-^"^ (6) 
/7(7,7)d7 = ^7(74,7), 

where /^(■) and Py(-) are the PDF and CDF of 7, respectively. 

Note that, according to the approximation model, the analytical FER can be calculated as 
an outage probability, and thus, the accuracy is mainly determined by the SNR threshold 
selection. 




B. Some Existing SNR Threshold Approaches 

In [|27ll . El Gamal and Hammons demonstrated the SNR threshold-based FER approxima- 
tion model for iteratively decoded systems employing turbo codes. And the optimal SNR 
threshold has been proved to coincide with the convergence threshold of the iterative turbo 
decoder. 

Recently, in [|28ll . [|29ll , Chatzigeorgiou extended this model to non-iterative coded and 
uncoded systems. To get the optimal SNR threshold, proper error criterion should be used. 
In [|28l . II29I . the minimum absolute error sum criterion is adopted to minimize the sum of 
absolute error 

7i = argmin|^ |P/(7) - P^(7t, 7) | d7 
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where 7 is the average SNR, and the SNR threshold can be calculated as 

However, the model on the basis of the minimum absolute error sum criterion might not 
be sufficiently accurate since it does not consider the fact that FER decreases more quickly 
at high SNR region in high diversity order systems. Hence, it can be improved. 



C. Model improvement and the SNR threshold result 

In this subsection, we propose an improved SNR threshold-based FER approximation 
model for the general diversity systems. The improvements are twofold. 

Firstly, taking into account the fact that the FER decreases quickly when SNR increases, 
minimum absolute relative error sum criterion, which minimizing the sum of absolute relative 
error, can be adopted 

P/(7)-F,(7i,7) 



7( = arg mm 



PAi) 



d7 



(8) 



As Eq. dH]) is difficult to solve, we instead use a suboptimal error criterion. Reasonably, 



suppose min < 



d7| < 00 and notice that the integration in Eq. ([8]) is from 
to 00, and then, the absolute relative error should approach zero when 7 — )• 00: 

^/(7)- ^7(7^,7) 



It 



arg lim 

7—^00 







(9) 



^ arg lim {P^(7) - ^,(7^, 7) = O} . 

7— !>00 

Otherwise, for a sufficiently big value T (0 < T < 00) and a small enough value 5 



(0 < 5 < 00), the absolute relative error can be greater than 5, i.e. 



P!{l)-Pi{lUl) 
Pf(l) 



> 



S, when 7 > T. Hence, the absolute relative error sum can't be minimized as it approaches 
infinity: 



Pf(l)-Fi{lt,l) 
Pjil) 



d7 > 5dq = 00. 

Secondly, the SNR threshold should include the factor of the diversity order of general 
diversity systems. For example, in a Rayleigh fading channel, the CDF of SNR is -F^(7, 7) = 
1 — e~^l^ , then lim^_>.oo ^7(7, 7) < 7"^7 and in a Nakagami-m fading channel, the CDF of 
SNR is given by: -F^(7, 7) 

the lower part incomplete Gamma function given by: r(m, x 

■mr(m)(7)™7 



^^"r'(my^'* ' where r(m) is the Gamma function, and r(m, x) is 

t™^^e~*dt, then using the 



fact that limt^o e^* < 1, we can get lim^^oo ^7(7, 7) < V{m)-^ Q{^)'^-^dit ^ - "^ 
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Suppose for some diversity system with a diversity order of d, the CDF of SNR can be 
approximated in a form Il29ll - [l3n : 

limF^(7,7)^G'(7)7', (10) 
7—^00 

where G{'j) is a constant related to 7. 

Then, combining Eq. ([5]), Eq. ^ and Eq. (flOl) . we can get the analytical SNR threshold 
for our improved model as 

7i,d ^ arg_lim |^ f,{j,^)P^{^)dj - ^(7)7' = o| 
= arg lim ( H G{^)dY~'Pf{l)d^ - G{^hf = 
= arg_lim |7/ = rf^ 7'^-ipf (7)d7| 



(11) 



{d ^^-'Pf{^)d^^ 



l/d 



In uncoded systems, Pf{'y) can be given in closed-form. For example, for linear modula- 
tion, Pf{'^) can be approximated as Il28]| 

Pf{^)^l-{l-Q{^))\ (12) 



/ 

where c is modulation constant (c = 2 for binary-phase-shift-keying (BPSK)), Q(-) is the Q 
function, L is the frame length. Substituting Eq. (fT2l) into Eq. (fTTI) . we can get SNR threshold 
7t_rf for uncoded systems: 



lt,d 



(^d Y''(l~{l-Q{^)f)A^^ . (13) 

In coded systems, we can calculate Eq. (fTTI) using numerical methods. We can first get 
P^(7), i.e. the instantaneous FER for the scheme over AWGN channel, using Monte Carlo 
methods. When the SNR values 7^, i = 1, 2, ■ ■ ■ , A^, are equally spaced with A7 and ordered, 
the following equivalent expression for discrete SNR values can be obtained 

/TV \ l/'i 

lt4^ [dY.lt 'Pf{ii)M\ ■ (14) 

Note that when c = 2 (BPSK), L = 1, and = 1, we can obtain 7t,i = j using Eq. (fT3l) 
[|32l . and P/ = 1 — ~ ^, which is equivalent to the well known average BER of BPSK 
over Rayleigh fading channel at high SNR: ^^ = 1(1-^^-)^^ [i26ll. 

If we consider a multiple-input multiple-output (MEMO) channel having N-r inputs and 
Nji outputs. The transmitter uses space-time block coding Il33ll . while the receiver coherently 
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combines the = NtNr independent fading paths. If 7 now corresponds to the instantaneous 
SNR at the output of the combiner, its probability distribution is given by [|26ll . [|29ll : 



{^/NTr{N-l)V 

where 7 is the average SNR per receive antenna. Using the fact that e^"^ < \ {x > 0), the 
CDF of 7 can be approximated as: 

1 



where G{^) = (^/^v^^pvjvT- So, we can obtain the SNR threshold 'jt^N using above method, 
and the approximated FER of the system for MIMO quasistatic fading channels can then be 
approximated as [|26l, |l29l 

Pf ^ F^(7t,7v) 



k\ 

k=0 

{^/Nt)^N\ 

So far, we have developed an improved SNR threshold-based FER approximation model 
for general diversity systems. In this model, FER is approximated as an outage probability 
as Eq. Q, the SNR threshold is given by Eq. (fTTl) and can be calculated using Eq. (fT3l) for 
uncoded systems and using Eq. (fT4l) for coded systems. 

IV. FER Analysis of the HRS Scheme 
In this section, we perform the FER analysis of the HRS scheme in cooperative systems. 

A. Analytical FER Analysis of the HRS Scheme 

According to the proposed FER approximation model, a frame error only occurs if SNR 
is below the SNR threshold ^t,d- Then, the average FER of the HRS scheme, denoted by Pf, 
can be expressed as 



(15) 
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where 'Jhrs is the SNR of the HRS scheme and d is the diversity order of the cooperative 
system with the HRS scheme. 

Hence, to get the FER for the HRS scheme, we merely need to derive the CDF of 'Jhrs- 
Since the channels are Rayleigh block fading, the instantaneous SNR is exponentially 
distributed. Then, the CDF of 'ju is given by 



K^ht,l) = 1 - e 



-Aii7t,i 



where Ah = I/7H is defined in Eq. dH). 

For simplicity, we introduce a vector variable z, where 



Z — [zi,Z2, - ■ ■ ,Zn],Zi 



0, if 7h > 7i,i, 

1, if 7h < 7t,i- 

In the HRS scheme, the i-th relay node is assigned to either Gdf or Qaf according to its 
SNR 7h. If 7ii > 7t,i (Zi = 0) then i E Gdf, and if 7^ < 7^,1 (Zi = 1) then i E Qaf- By 
using the fact that for any value x, a;° = 1, then the probability of Zi, denoted by Pr{zi), 
can be written as 

Pr{zi) = [Pr{zi = 0)]'-'^[Pr{zi = 1)Y^ 

= [Pr{ju>lt,i)Y-''[Prhu<lt,iT 

Then, the CDF of 'Jhrs is derived in Appendix |A] as 

FjHRsht,d) ~ 2^ Ci— r — r ^ ^^^^ 

where b = {6„ ^ = l, . . . , n} g Z„, 6,, G {0, 1}, Ci = (-1)^7=1 and C2 = Er=i H^,>^ii + 

Substituting Eq. (fT6l) into Eq. (fTSl) . we can obtain the FER for the HRS scheme 

^ ^ C:,(l-e-Ao7M)-Ao(l-e-^^^M) 

5. FEi? Simplification at High SNR 

In this subsection, we derive the simplified FER for the HRS scheme at high SNR. At 
high SNR, To, 'fu and 7^j are large, then according to Eq. (HI): Aq = I/tq, Aij = l/'fu and 
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= 1/72* ^re small. Using the approximation that 1 — exp(— x) ~ x, when x — )■ 0, we can 
get the simplified CDF of 'Jhrs at high SNR is given by 

\ n+l n 

F,,ns{lt4) ^ ^fVn( — + (18) 

The proof of Eq. (fTSi) is given at appendix |Bj 

Substituting Eq. (fTSl) into Eq. (fT5l) . we can obtain the simplified FER of the HRS scheme 
at high SNR 

p/-^^AonAji+^^y (19) 

n + l J-JL \^ 1t,d^2iJ 

2=1 ' ^ 

where Aq IliLi (-^li + ^21) (7)^'-""*'^^ which indicates that the diversity order of n + 1 can 
be achieved, d = n + l, 'jt^^ = 'jt,n+i and see also Eq. (I26l) for the last step. 

The simplified FER is more intuitive and simple compared to the FER expression in Eq. 
(fTTl) . Note that both Eq. (flTl) and Eq. (fT9l ) are analytical approximate average FER expressions 
for the HRS scheme in cooperative networks. The calculation methods of parameters 74 1 and 
7t,rf are presented in Section |inl 

C. Extending to other Relaying Schemes 

The FER results for the HRS scheme in Eq. (fTTl) and Eq. (fT9l ) can be easily extended to the 
cooperative systems using the AF with relay selection (AF-RS) scheme [|34l and the perfect 
DF with relay selection (PDF-RS) scheme [[35l. In the AF-RS scheme, a single relay with the 
maximum source-relay-destination SNR will be selected to forward signals from the source 
to the destination using AF [34J. In the PDF-RS scheme, it assumes that all relay nodes can 
correctly decode the received signals and a single relay with the maximum relay-destination 
SNR will be selected to forward signals from the source to the destination. For simplicity, 
we only present the simplified results at high SNR [35|. 

For PDF-RS, the transmissions from the source to any relay are always correct, it means 
that 7ij = 00 and Zi = 0. Thus, the average FER of PDF-RS, denoted by p^^^-^^^ is given 
by 

n + l 

i=i 

For AF-RS, the SNR of the z-th Rayleigh relay link (from the source to the z-th Rayleigh 
relay and to the destination) can be approximated as "combined" Rayleigh fading with 
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Fig. 2. FER comparison of proposed model and the model of P281 for case 0: general MIMO channels with A'^r = 1, 

N ^ Nr ^ 1,2,4, uncoded. 

parameter of Xf^ ^ Xu + X2i at high SNR Q. Using this approximation, AF-RS can be 
viewed as a special case of PDF-RS with X'^i = Xf^. The average FER of AF-RS, denoted 
by p^^-^^ ^ is given by 



P 



i=i ^ ' 



Defining the gain of the HRS scheme over the AF-RS scheme as the ratio of FER, we can 
get 

pAF-RS n 1 _|_ All 

G HRS = ^ TT -^(dB). 

AF-RS p. ii]^ I TmAu^ ^ 

It indicates that the gain improves when ^ increases, and when the number of relay nodes 



increases, the gain increases quickly. 



V. Simulation results 

In this section, we compare the proposed approximate FER expressions with the exact 
FER obtained by Monte-Carlo simulations. Unless specifically mentioned, the simulations 



'This is a lower bound for AF 1361 . 1371 . 
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Fig. 3. Average FER of the HRS scheme for case 1: Jlo = ^lu — Q,2i ~ 1, uncoded. 
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Fig. 4. Theoretical and simulated FER of the HRS scheme for case 4 and 5: Jlo = ^^u ~ = 1, coded. 
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A n=4, case 7, simulated. 
— n=4, case 7, theoretical. 
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Fig. 5. Theoretical and simulated FER of the HRS scheme for case 6 and 7: fio = ^^u ~ = 1, coded. 
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Fig. 6. Average FER of the HRS, AF-RS and PDF-RS schemes for case 4 with n = 4: fio = Qu = fl2i = 1, coded. 
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TABLE I 
SCENARIOS 





System 
Type 


Frame 
Size 


uncoded/ 
coded 


Code 
Rate 


Number of 
Nodes 


oiNKS 


case 


MIMO 


100 


uncoded 


- 


Nt=1, N=Nfl=l,2,4 


fl = l 


case 1 


HRS 


100 


uncoded 




n=l, 2, 4 


Qo — Qli = ^l2i ~ 1 


case 2 


HRS 


100 


uncoded 




n=l, 2, 4 


rif) = 1, iiii = 16, r22i = 1 


case 3 


HRS 


100 


uncoded 




n= 1, 2, 4 


fio = 1, nu = 1/16, !r22« = 1 


case 4 


HRS 


100 


coded 


1/2 


n= 1, 2, 4 


Qo — Qli = Q2i ~ 1 


case 5 


HRS 


200 


coded 


1/2 


n= 1, 2, 4 


fio = Qli = f22i = 1 


case 6 


HRS 


100 


coded 


2/3 


n= 1, 2, 4 


flo = fill = S^2i ~ 1 


case 7 


HRS 


200 


coded 


2/3 


n= 1, 2, 4 


fio = ^^li = Q2i ~ 1 



are performed for a BPSK modulation and a frame size of 100 or 200 symbols over block 
Rayleigh fading channels. We consider two basic cases: (1) the uncoded case, where none 
channel code is used; and (2) the coded case, where a systematic convolutional code with a 
code rate of 1/2 and the generator matrix of (5, 7)8 or a code rate of 2/3 and the generator 
matrix of (23, 35, 0; 0, 5, 13)8 is used. 

It is practically infeasible to verify the accuracy of analytical results for all the possible 
scenarios, as a prohibitively large number of combinations can be generated by varying relay 
number n and average SNR parameters Therefore, the scenarios in Table |I] are chosen 
such that a range of diverse cases are covered. For instance, cases 1, 2, and 3 are uncoded 
and case 4, 5, 6 and 7 are coded. Cases 1, 4, 5, and 6 are symmetric situations, while cases 
2 and 3 are dissimilar situations. In all cases, 1, 2 and 4 relay nodes are considered. We also 
choose case for comparing our proposed FER approximation model to the results of [28 1. 
And in all cases, = 'jt,n+i- The SNR threshold values are obtained using the methods 
developed in Section [nil and are shown in Table [III 

Fig. [2] shows the analytical FER curves by using the proposed model and the model of 
||28]| for case 0. For case 0, the SNR threshold was found to be 4.6 dB based on Eq. ([7]) 
(the model of [|2l), and was found to be 5.10 dB, 5.36 dB and 5.89 dB for = 1,2,4, 
respectively, by our proposed model based on Eq. (fT2l) . As shown in the figure, our proposed 
model converge with the simulated FER quickly as the SNR increases while the relative error 
using [|28ll cannot be ignored even at high SNR. The FER results based on the model in [|28l 
become less accurate when the diversity order increaeses, but our proposed model is still 
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TABLE 11 
SNR Threshold 





System Type 


uncoded/coded 


Number of Nodes 


Diversity Order 


SNR Threshold 








Nt=1, N=Nfl=l 


1 


5.10 dB 


case 


MIMO 


uncoded 


Nt=1, N=Nfl=2 


2 


5.36 dB 








Nt=1, N=Nfl=4 


4 


5.89 dB 








n=l 


2 


5.36 dB 


case 1,2,3 


HRS 


uncoded 


n=2 


3 


5.62 dB 








n=4 


5 


6.16 dB 








n=l 


2 


-1.48 dB 


case 4 


HRS 


coded 


n=2 


3 


-1.26 dB 








n=4 


5 


-0.81 dB 








n=l 


2 


-0.85 dB 


case 5 


HRS 


coded 


n=2 


3 


-0.68 dB 








n=4 


5 


-0.33 dB 








n=l 


2 


0.06 dB 


case 6 


HRS 


coded 


n=2 


3 


0.20 dB 








n=4 


5 


0.47 dB 








n=l 


2 


0.60 dB 


case 7 


HRS 


coded 


n=2 


3 


0.70 dB 








n=4 


5 


0.90 dB 



accurate. While not shown here, similar trends can be observed for other scenarios. 

Fig. [3] shows the proposed approximate FER results and simulated results of the HRS 
scheme for case 1 . From the figure we can see that the analytical theoretical curves obtained 
by Eq. (flTl) converge well with the simulated results, and the simplified results computed 
by Eq. (fT9l) converge to the simulated results at high SNRs. Figs |4] and [5] show the different 
encoders and other block lengths for coded scheme and confirm the good performance of 
the proposed model for different parameters. Case 2 and case 3 are also verified through 
simulation but omitted here for brevity. 

Fig. [6] shows our proposed simplified FER of the AF-RS and PDF-RS for case 4 with 4 
relay nodes. From the figures we can see that the simplified results converge to the simulated 
results at high SNRs for both AF-RS and PDF-RS. And the HRS scheme has considerable 
performance gain over AF-RS. 
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VI. Conclusions 



In this paper, we have analyzed the average FER of the HRS scheme in general cooperative 
wireless networks with and without applying channel coding at transmitting nodes. We 
proposed an improved SNR threshold-based FER approximation model. We then apply this 
model to HRS system and derived the analytical approximate average FER expression and the 
simplified asymptotic FER expression at high SNRs for the HRS scheme. Simulation results 
match well with the theoretical analysis, which validates our derived FER expressions. 



According to Eq. ([3]), 'Jhrs = 7o + 7m- In the following, we first derive the PDF of 70 
and 7m and then get the PDF of 'Jhrs as the convolution of the PDF of 70 and 7^. 
As 7o, 7ii and are exponentially distributed, their PDF and CDF are given as 



Appendix A 



The PDF and CDF of -/hbs- The Proof of Eq. ([16]) 



,ke {0,li,2i}, 



and 



^7.(7) = l-e 



,ke {0,li,2i}, 



respectively. 

According to Eq. ([T]) and Eq. 7^ = maxi<j<„{7j} and 




To get the distribution of 7^, we first derive the distribution of %. 
If Zi = 0, 7j = 72j. Then, the conditional CDF of 7j given = is 



F^^ (7 1 = 0) = Pr (72i < 7) = 1-e 



(20) 



If Zi = 1, using the approximation 
7i given Zi = lis [IH. 



x+y+l 



xy_ 



min{x,y} [|36l . then the conditional CDF of 



^7»U,=i(7) = Pr{min{ju\,^=i,-f2i} < 7) 

= 1-(1-F,,,|,^^,(7))(1-F,,,(7)). 



(21) 
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The conditional PDF and CDF of 'ju\zi=i are 

r I X f-iMP'rinii = l,Zi = 1) 
J7iiU,=il7j = 



Pr{z, = 1) 

0, if 7 > 7t,i, 

(22) 

1, if 7 > 7t,i' 

where Pf^^^{l) is the FER of 7ij at AWGN channel and Pf.^^^{l) ~ 1 if 7 < 7t,i and 

^/:7u(7)~0if7>7t,i. 

Combining Eq. (|2TI) and Eq. ((22l) . the corresponding conditional CDF of 7^ given Zi = 1 
can be written as 

f 1 - ^-'^'^-f V-A2.7 if 7 < 7, 1 
F^,|.,=i(7) ~ < (23) 
[1, if7>7t,i, 

where at the second step, we extend the range of 7 from 7^ 1 to 7^^ using an exponential 

function to approximate F^^i^— 1(7). The simulation results in Section IVl indicate that this 

approximation is accurate. 

Combining Eq. (|20l) and Eq. (|23]) and noting that Pr{zi = 1) = 1 - e"^"''".! and Pr{zi = 

0) = e^'^"'''' ^ the CDF of 7j can be rewritten as 

F,^(7) = Pr{zi = l)F,,(7k. = 1) + Pr{zi = 0)F,Xl\z^ = 0) 

1 _ e-(^i>+^207^ if7<7t,i, 

1 - e-(^i'^'-i+^2.7)^ if ^ > ^^^^^ (24) 

1-e ^t.d , if 7 < 7i,d, 

1 _ e-(Ai,7t.i+A2,7)^ if7>7t,d, 

Noting that 74 1 ^ 74 ^, we can use a exponential function to approximate the CDF of 7^ 
when 7 < 7t^rf. In order to keep the value Fj^{^t,d) unchanged, we have the last step. Then 
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the CDF of 7m (when 7 < 7^^^) is given by [|38l 

^7m(7) = Pr{m&]i{li} < 7) 



l<i<n 



1-e (25) 



i=l 



with the help of that 

n n n n n 

i=l 1=1 1=2 i=l b£Zn 4=1 

Where b = {6„ ^ = 1, . . . , n} e Z„,, 6, e {0, 1}, C^ = (-1) and C2 = ELi H^,>^ii + 
The PDF of 7^ is given by 

be2„,b^0 

As -fHRs = 7o + 7m, the PDF of 'Jhrs when 7 < 7^,^ can be expressed as 

ri 

flHRS (7) = / /70 (7 - t) {t)dt 







- rAoe-'''(^-*)5^CiC2e-'^^*dt 
-Aoe-^°^ J] C1C2 / e-^^'-^°'^'dt 



(27) 



1 - e-('^2-'^o)'>' 



C-2 — An 

be2„,b^o ^ " 

g-Ao7 _ g-C27 



5^ A0-C2 

be^n.b^o 



Finally, the CDF of Jhrs is given by 

rit.d 





5Z ^1^°^^ A0-C2 (28) 
" be^„,b^o 



^ ^ W-e- ^"^^^'^) - Ao(l - e-'^^'^'.'*) 



An — C' 

be2„,b7^0 ^ 
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Appendix B 

Simplified CDF of ^hrs at high SNR: The Proof of Eq. (fTSl) 

Using Eq. (|27]) and the approximation that e"^ ^ 1 — x, the PDF of 'Jhrs at high SNR 
can be approximated as 

(1 - Ao7) - (1 - C27) 



Aq — C2 



(29) 



^7,.(7)^7"n(— + (30) 



heZn 

= -Aq ^ C1C27. 

Applying the approximation that 1 — x to Eq. (|24l) and Eq. (|25l) . we can get 

n 

.7t,i 

i=l 

and 

F,^il) ^ E "^i^l - ^27) = E - E '^i^27- (31) 

heZn hez„ hez„ 

As is CDF, F^,„(0) = 0, we can get 

F,^{0)=J2Ci = 0. (32) 

h£Z„ 

Combining Eq. (I30l) . Eq. (|3TI) . and Eq. (I32l) . we can get 

n 

- E C1C27 ~ 7" n^— + ^2^- (33) 

heZn 'i = l 



Substituting Eq. (1331) into Eq. (I291) . we can get the approximated PDF of 'Jhrs as 

n 

.7i,i 

And the CDF of 'Jhrs can be expressed 



7i,rf 



nt,d An7"j"^ " 'y 

Jo a -h i ^t,d 
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